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PURPOSE 


The  purpose  of  the  project  is  to  increase  understanding  of  the  roles  of 
cloud  conversion,  accretion,  evaporation,  and  entrainment  processes  in  shaping 
the  distributions  of  water  vapor,  cloud,  and  precipitation  associated  with  tropical 
circulations. 


ABSTRACT 

A  system  of  conservation  equations  embodying  conversion  of  cloud  to 
precipitation  and  cloud  collection  by  precipitation  are  examined.  The  steady 
precipitation  rate  at  the  base  of  a  model  updraft  column  is  shown  to  be  a 


maximum  when  cloud  conversion  and  collection  are  complete  above  the  level 

of  nondivergence  and  zero  below  that  level.  The  decrease  of  saturation  vapor 
"  ■  /  .. 
density  with  height,  characteristic  of  the  atmosphere,  implies  that  for  any 

cloud-water  content,  there  is  a  level  (compensation  level)  above  which  the 

cloud  density  decreases  with  ascending  motion,  and  conversely.  For  a 

particular  choice  of  the  magnitude  of  cloud  collection  1  the  precipitation  rate 

at  the  ground  and  its  time  of  onset  is  rather  insensitive  to  the  rate  of  cloud 

conversion.  The  efficiency  of  precipitation  production) at  the  earth's  surface 

by  transient  disturbances  in  a  saturated  model  atmosphere  decreases  with 

increasing  speed  of  the  updraft.  Paths  for  further  study  of  the  role  of  cloud 

conversion  and  cloud  collection  processes  for  the  initiation  of  precipitation 

are  indicated.  The  difference  between  the  horizontal  speed  of  raindrops  and 

the  horizontal  wind  speed  is  almost  always  less  than  1%  of  representative 

horizontal  winds,. 


PUBLICATIONS,  LECTURES,  REPORTS,  AND  CONFERENCES 


A  paper  entitled  "Relationships  Between  Tropical  Precipitation  and 
Kinematic  Cloud  Models"  was  submitted  to  and  accepted  by  the  Program 
Committee  for  the  43rd  Annual  Meeting  of  the  American  Meteorological 
Society  (New  York  City).  This  paper  is  scheduled  for  presentation  on 
January  21,  1963. 

On  September  4,  Drs,  Kessler  and  Newburg  and  Mr.  Feteris  visited  the 
U.S.  Army  Electronics  Research  and  Development  Laboratory  and  discussed 
various  matters  pertaining  to  the  work  in  progress  and  future  plans. 


1.0  FACTUAL  DATA 

1,1  A  Model  of  the  Bulk  Effects  of  reoipitation-Cloud  Interactions 

The  derivation  of  continuity  equations  for  cloud  and  precipitation  and 

exposition  of  methods  for  solving  them  are  outlined  in  Sections  1.4ff  of  the 

first  report  of  this  contract,  and  background  details  are  in  the  literature 

[7,  8,  9].  In  their  general  forms,  the  equations  considered  are 

8M  9M  9M  9M  9  „  91np 

—  =  -u— -V  —  -  W  —  -  —  MV  +  Mw 

8t  9x  8y  9z  Bz  8z 

+  cloud  conversion  +  cloud  collection  -  evaporation 


8m  _  _  9m 
9t  U  9x 


9m  9m  81np 

v  — —  -  w  r— 1  +  wG  +  mw 
9y  9z  9z 


-  cloud  conversion  -  cioud  collection  +  evaporation. 

Equations  (1)  and  (2)  incorporate  the  following  continuity  equation  for  air: 

+  +  if  +  W  =  °*  (3) 

In  Eq.  (l)r  M  (always  2:0)  is  the  density  of  precipitation  in  grams  per  cubic 
meter  and  refers  to  that  portion  of  the  condensate  which  has  an  appreciable  fall 
speed  relative  to  the  air.  M  is  assumed  to  be  distributed  over  an  exponential 
population  (Marshall-Palmer  distribution  [12])  of  particle  sizes.  V  <  0  is  the 
terminal  fall  speed  of  the  precipitation  particle  whose  diameter  DQ  divides  the 
exponential  distribution  into  parts  of  equal  water  content.  The  components  of 
air  motion  in  the  x-  and  y-  (horizontal)  and  z-  (vertical)  directions  are  given  by 
„  u,  v,  and  w  respectively;  p  is  the  density  of  air;  and  t  is  time. 

In  Eq.  (2),  m  is  the  cloud  density  plus  the  vapor  density  minus  the  satura¬ 
tion  vapor  density,  and  may  be  positive,  negative,  or  zfiro.  When  m  is  positive, 
it  represents  the  density  of  cloud,  defined  as  condensate  whose  terminal  fall 
speed  is  zero.  When  m  is  negative,  it  represents  the  amount  by  which  the  actual 
.vapor  density  is  less  than  saturation.  G,  a  function  of  height,  is  called  the 

generating  function:  G  =  -  p  (dQ  /dz),  where  Q  is  the  Saturation  mixing  ratio 

s  a 


i 


of  water  vapor  In  air.  The  cloud-conversion,  cloud-accretion,  and  evaporation 
terms  are  discussed  in  Report  No.  1;  they  define,  respectively,  the  spontaneous 
conversion  of  cloud  to  precipitation,  the  collection  of  cloud  particles  by  precip¬ 
itation,  and  the  evaporation  of  precipitation.  Because  the  model  requires  immediate 
evaporation  of  cloud  in  subsaturated  air,  evaporation  of  precipitation  does  not 
occur  at  the  same  time  and  place  as  cloud  conversion  and  accretion. 

Equations  (1)  and  (2)  relate  local  time  changes  of  precipitation,  cloud,  and 
vapor,  to  transport,  generation,  and  exchange  processes.  The  first  three  terms 
on  the  right  of  both  equations  measure  effects  of  advection  by  the  wind.  The 
term  9(Mv)/3z  in  Eq.  (1)  measures  the  vertical  divergence  of  the  mass  flux  of 
precipitation.  This  term  does  not  appear  in  Eq.  (2)  because  V  =  0  for  cloud. 

The  terms  in  w[(81np)/8z]  account  for  the  compressibility  of  a  horizontally 
uniform  atmosphere.*  Equations  (1)  and  (2)  can  be  solved  for  time-dependent 
fields  of  precipitation,  cloud,  and  vapor  when  initial  and  boundary  conditions  with 
the  wind  field,  the  generating  function,  and  the  cloud  physics  terms  are  specified. 


1.2  Idealized  Steady-state  Precipitation  Rates  from  a  Saturated  Horizontally 
Uniform  Updraft  Column 

At  a  shower  core  or  near  the  center  of  a  widespread  precipitation  abea, 
the  horizontal  advection  terms  in  Eqs.  (1)  and  (2)  vanish  and  the  sum  of  the 
equations  is  “ 


3  (m+M)  _ 

8t  =  "  w 


8(M+m) 

3z 


~  MV  +  wG  +  (m  +  M) 

0  Z  0Z 


(4) 


Consider  the  vertical  integral  of  Eq.  (4)  in  the  steady  case: 


v 


0Z 


*When  M,  m,  and  G  are  in  mixing  ratio  units,  Eqs.  (1)  and  (2)  are1  the  same 

except  that  Mw[<81np)/8z]  in  Eq.  (1)  becomes  -MV[(81np)/8z] ,  mw[(81np)/8z] 

is  completely  dropped  from  Eq.  (2),  and  G  =  -dQ  /dz. 
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The  limits  of  integration  are  from  the  base,  z  =  0,  to  the  top,  z  =  H,  of  the 
updraft  column.  The  second  integral  on  the  right  is(MV)  -  (MV)  ;  since 

Z— Jtl  Z™  u 

M  *  0  at  the  top  of  the  updraft  column,  this  is  just  the  precipitation  rate  RQ  at 
the  base  of  the  updraft.  can  be  taken  as  the  precipitation  rate  at  the  ground 
without  loss  of  generality.  The  third  integral  on  the  right  is  the  condensation 
rate  in  a  vertical  column  of  unit  cross  section. 

The  first  integral  is  studied  by  integration  by  parts.  Note  that 

"  w  =  “  9-7  fw  <m  +  M)]  +  (m  +  M)  fj.'  .  (6) 

Tile  integral  from  z  =  0  to  z  -  H  of  the  first  term  on  the  right  of  Eq.  (6) 
vanishes  because  w  is  zero  at  the  base  and  at  the  top  of  the  updraft  column.  The 

divergence  of  the  vertical  wind  in  the  second  term  on  the  right  can  be  replaced 

»  . 

by  its  equivalent  in  the  equation  of  continuity  for  air,  Eq.  (3).  Substitution  from 

•> 

Eq.  (3)  into  Eqp.  (6)  and  (5),  with  the  other  simplifications  noted  above,  gives 


Equation  (7)  is  a  reminder  that  the  horizontal  divergence  of  the  wind  is 
implicit  even  in  the  one-dimensional  forms  of  Eqs.  (1)  and  (2).  The  last  term 
in  Eq.  (7)  may  be  positive,  negative,  or  zero,  and  is  largely  dependent  on  the 
magnitude  and  distribution  of  the  ratio  of  precipitation  fall  speed  to  the  updraft. 
Because  horizontal  convergence  characterizes  the  lower  part  of  an  updraft 
column,  and  horizontal  divergence  the  upper  part,  precipitation  rates  in  excess 
of  the  condensation  rate  are  favored  by  relatively  large  values  of  m  and  M  below 
the  level  of  no  horizontal  divergence."  Equation  (7)  is  the  same  whether  the  atmos¬ 
phere  is  assumed  to  be  compressible  or  incompressible.  However/ as  shown  else¬ 
where  for  the  cloud-free  model  [7,  Fig.  1]  and  for  m  in  Section  1.3  of  this  report, 

the  vertical  distributions  of  M  and  m  vary  somewhat  with  the  considerations  of 

\ 

compressibility,  and  so  must  the  effect  of  the  divergence  term  in  Eq.  (7). 
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It  is  possible  to  treat  certain  simple  distributions ‘of  cloud  and  precipita¬ 
tion  by  exact  methods  paralleling  those  described  in  the  references.  These 
cases  are  interesting  as  limiting  forms  associated  with  limiting  values  of  the 
precipitation  rate.  They  indicate  how  extreme  modification  of  cloud-precipita¬ 
tion  interactions  could  affect  the  spatial  distributions  of  cloud  and  precipitation 
and  the  areal  distribution  of  the  precipitation  rate  at  the  ground. 

Case  1.  The  cloud  conversion  term  is  identically  zero.  When  the  cloud 
conversion  term  is  identically  zero,  M  is  zero  throughout  the  depth  of 
the  updraft  column,  and,  therefore,  so  must  be  the  precipitation  rate  at  the 
ground.  In  this  case,  manipulations  like  those  used’ to  derive  Eq.  (7) 
readily  show  that  the  terms  on  the  right  of  Eq.  (7)  are  equal  and  opposite  in 
sign.  In  this  case,  condensate  existing  in  the  updraft  column  is  removed 
by  divergence  aloft  as  rapidly  as  it  is  formed.  The  cloud  either  evaporates 
elsewhere  in  downdrafts  or,  after  changing  to  precipitation,  falls  toward 

;;  ..  ’  .....I/ 

the  ground  at  some  place  removed  from  the  updraft  core.  The  steady 
vertical  distribution  of  cloud  in  an  incompressible  atmosphere  is  given 

by  .  ..  ..  -  .  ..... 

a  m(z)  =  Jq  G  dz .  "  (8) 

Figure  1,  top,  shews  the  steady  vertical  distribution  of  cloud  for  case  1 
for  constant  and  linear  distributions  of  G  in  an  incompressible  atmos¬ 
phere.  As  Eq.  (8)  shows,  the  cloud  distributions  are  independent  of  the 
updraft,  distribution. : 

Case  2.  The  cloud  conversion  term  is  zero  above  the  level  of  nondivergence 

H  ' 

below  this  level  cloud  conversion  and  collection  by  precipitation  proceed  so 

"  •  '  .  ii  " 

rapidly  that  no,  cloud  coexists  with  precipitation.  ,  In  this  case,  the  precip¬ 
itation  content  is  zero  at  the  level  of  nondivergence  and  increases  in  the 
downward  direction  from  there  to  the  ground.  Above  the  level  of  non¬ 
divergence,  all  condensation  remains  as  cloud  and  is  borne  away  by 
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horizontal  divergence.  The  distribution  of  cloud  above  the  level  of 
nondivergence  is  given  by  an  equation  like  (8).  The  distribution  of 
precipitation  below  the  level  of  nondivergence  is  given  by 

“W  =  dz-  <9> 

H/2  • 

Case  2  distributions  of  cloud  and  precipitation  are  shown  in  Fig.  1, 

center,  for  two  ratios  of  precipitation  fall  speed  to  maximum  updrafts 

and  for  two  distributions  of  G  in  an  incompressible  atmosphere.  The 

2 

distributions  of  M  are  based  on  updrafts  w  =  [(4w  )/(H)](z  -  z  /H) 

max 

in  this  and  following  cases. 

Case  3.  The  cloud  conversion  and  collection  terms  are  so  large  that 
all  condensate  appears  immediately  as  precipitation  of  fall  speed  V. 

This  is  the  cloudless  case  discussed  in  the  references.  For  constant  V, 
the  steady  surface  precipitation  rate  at  updraft  cores  always  exceeds 
the  vertically  integrated  condensation  rate  in  a  vertical  column.  As 
the  ratio  of  precipitation  fall  speeds  to  updrafts  declines  toward  unity, 
the  core  precipitation  rate  increases  at  the  expense  of  the  rates  away 
from  the  core.  On  the  other  hand,  when  fall  speeds  are  much  larger 
than  maximum  updrafts,  the  core  precipitation  rate  is  only  slightly 
larger  than  the  vertically  integrated  condensation.  Distributions  of 
precipitation  in  this  case  are  shown  in  Fig.  1,  bottom,  for  two  ratios 
of  precipitation  fall  speed  to,  maximum  updrafts  and  for  two  distribu¬ 
tions  of  G  in  an  incompressible  atmosphere. 

Cases  4  and  $.  Above  the  level  of  nondlyergence,  cloud  conversion  and 
collection  proceed  so  rapidly  that  cloud  is  nonexistent;  below  the  level 
of  nondivergence,  conversion  and  collection  are  zero,  and  at  each  height 
the  cloud  amount  equals  the  amount  of  condensation  in  air  lifted  from 
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Cloud  and  precipitation  distributions  in  Idealized  one-dimensional  cases* 


Table  4b 
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the  surface.  These  conditions  with  V  =  constant  comprise  case  4; 
the  same  with  a  five  fold  increase  in  falling  speed  at  the  level  of 
nondivergence  is  ease  5,  The  V-discontinuity  in  case  5  resembles 
that  accompanying  the  melting  of  snow.  This  model  distribution 
is  the  condition  of  maximum  precipitation  rate  at  the  updraft  base  “ 

for  specified  fall  speed  and  maximum  updraft.  Above  the  level  of 
nondivergence,  no  cloud  exists,  so  none  is  lost  by  horizontal  diver-  - 
gence;  below  the  level  of  nondivergence,  cloud  condensed  over-  a 
wide  area  is  swept  toward  the  updraft  center,  where  it  is  collected 
by  precipitation  at  the  level  of  nondivergence  and  borne  to  the 
ground.  Since  precipitation  cannot  exceed  condensation  over  a  great 
region,  excessive  precipitation  at  the  updraft  center  is  accompanied 
by  a  balancing  deficit  away  from  the  center.  The  one-dimensional 
cases  Imply  the  major  features  of  two-  and  three-dimensional  models 
required  to  satisfy  conservation  principles. 

The  computation  of  M-profiles  in  cases  4  and  5  involves  some  special 
considerations.  From  z  =  H  to  the  level  of  nondivergence,  M  is  given  by  an 
equation  like  (9).  From  z  =,  0  to  the  level  of  nondivergence,  m  is  given  by 
Eq.  (8).  At  the  level  of  nondivergence,  precipitation  M  falling  at  velocity  V  +  w 
collects  all  of  the  cloud  m  rising  at  speed  w.  Because  conditions  are  steady, 
the  value  of  M  below  the  level  of  nondivergence  is  defined  by  equating  the  trans¬ 
port  of  M  +  m  immediately  below  the  level  of  nondivergence  to  the  transport  Of 
M  above  that  level;  i.e.,  _ •  . 

<mw)below  *  ‘“(V  +  w,lbelow  '  +  w>l  above'  W 

The  updraft  w  has  the  same  value  immediately  above  and  below  the  level  of  non¬ 
divergence. 

Equation  (10)  has  been  applied  to  determination  of  the  discontinuous  M-profiles 
shoWh  in  Fig.  2,  where  the  upper  diagram  pertains  to  constant  G,  cases  4a  and  5a, 
and  the  lower  to  a  linear  distribution  of  G,f  cases  4b  and  5b.  Fdur  profiles  arc 


Fig.  3.  Schematic  illustration  of  .Cases  3  and  4,  discussed  in  the  text  and  listed  in 
Tables  3a,  b  and  4a,  b.  Horizontal  divergence  in  the  high  atmosphere  accompanies 
rising  air  motion  and  spreads  precipitation  packets  horizontally  as  they  descend  in  the 
shaded  area  marked  D.  The  packets  are  contracted  again  by  horizontal  convergence 
in  the  lower  atmosphere.  Precipitation  at  the  ground  in  the  section  marked  B  is  due 
to  condensation  in  a  volume  larger  by  Hie  extent  of  horizontal  shading  than  the  unshaded 
cross  section  vertically  overhead.  When  condensate  augments  precipitation  at  the 
level  of  nondivergence  but  not  below  that  height,  the  region  contributing  to  precipitation 
at  B  is  further  increased  by  the  area  of  vertical  shading  marked  A  in  the  figure.  Cloud 
in  area  A  is  svtept  inward  and  upward  by  horizontal  convergence  and  ascends  to  the 
level  of  nondlvergenoe,  where  it  is  collected  by  precipitation  and  ultimately  deposited 
at  the  ground  in  Section  B. 


( 
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given  in  each  diagram:  two  pertain  to  the  cases  in  which  V  has  constant  values  . 
of  -1  m/sec  and  -0.67  m/sec;  the  other  two  show  the  M-distribution  that  exists 

when  each  of  these  fall  speeds  increases  by  a  factor  of  five  at  the  level  of  non- 

•  |i 

divergence.  Such  a  discontinuous  change  of  V  approximates:  the  fall-speed 

changes  that  commonly  occur  at  the  melting  level. 

The  profiles  shown  in  Figs.  1  and  2  can  be  scaled  to  any  distribution  of 

atmospheric  parameters  that  differ  by  constant  factors  from  those  used  here. 

For  new  scale  height^  new  generating  functioiy</f  new  vertical  air  apeed^ty^J 

and  new  fall  speed  'If'  V/w  the  new  solutions  denoted  by/^are  given 
To  ■  /  'max  max  C  - 

by  wP/fin. 


Principal  features  of  the  profiles  in  Figs.  1  and  2  are  listed  in  Tables  1 
through  5,  which  are  numbered  to  correspond  to  case  numbers  of  the  figures. 
Note  that  Table  1  lists  values  of  w  and  /wQ  da,  which  are  the  same  for  all  the 
other  tables.  Note  that  the  surface-precipitation  rate  exceeds  the  vertical 
Integral  of  the  condensation //rate  in  practically  all  cases,  and  Is  more  than  five 
times  the  condensation  term  in  the  slow-fall-speed  example  of  case  4b.  Cases 
3  and  4  considerably  exaggerate  the  probable  magnitude  of  the  precipitation 
excess  that  characterizes  centers  of  natural  precipitation,  since  fall  speeds  In 
nature  are  usually  much  larger  than  updrafts,  The  mechanism  responsible  for 
the  very  high  precipitation  rates  of  case  4  is  schematically  Illustrated  In  Fig,  3. 
Table  6  shows  how  the  mechanism's  effects  are  minimized  when  the  fall  speed 
Increases  during  descent.  This  Is  discussed  in  some  detail  in  [7] .. 

1.3  Steady-state  Vertical  Distribution  of  Cloud  "in  a  Precipitation-free 
Compressible  Atmosphere 

The  maximum  possible  cloud  density  at  various  atmospheric  levels  is  a 
parameter  of  some  importance  to  considerations  of  precipitation-cloud  Inter¬ 
actions.  The  steady  distribution  of  cloud  established  in  an  air  column  after 
ascent  has  continued  for  a  very  long  time  without  fallout  of  condensation  products 

-  .  n  ii  " 

■  ;|t  -•  ■  \  -  '.  •  -  -  .. 
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also  represents  at  each  height  the  water  content  associated  with  ascent  from 
the  condensation  level  with  retention  of  all  cloud.  The  steady  profile  for  a 
given  condensation  level  thu^  represents  the  maximum  attainable  cloud  content 
at  each  height  for  that  condensation  level. 

In  the  steady  state  without  cloud  conversion  of  horizontal  advection, 

Eq.  (2)  . becomes 

91nc 


w 


0- 

am 


mw 


or 


dz  9z 

8m  31np 


-  wG  =  0 


m 


=  G 


(U) 


(12) 


8z  8z 

Sipce  w  does  not  appear  in  Eq.  (12)  the  solution,  as  expected,  is  independent  of  the 

shape  of  the  updraft  distribution.  We  assume  G  =  A  +  Bz,  and  take  alnp  =  K, 

'  ■  -  0Z 

a  constant.  These  are  reasonable  approximations  to  the  real  atmosphere.  Then 


the  solution  of  £q.  (12)  Is 


m  =  eK(Z  “  Z0)[m()  (A  +  BzQ  +  |)  ]  -  £  (A  +  Bz  +  j|), 


(13) 


In  the  tropical  atmosphere  A  *  3  X  10  3  gm/m4,  B  *  -3  X  10  7  gm/m5,  and 
-4  -1 

K  ®  -10  m  ;  substitution  of  these  values  in  Eq.  (13)  yields: 


m 


e10‘4<Z0 


Z)  [mQ  -  80  +  3  X  10  3  zfl]  +  [60  -  3  X  10~3z] . 

(14): 


-3 


The  units  of  m  are  gm  meter  when  z  is  given  in  meters. 


9  m 


The  location  of  maxima  of  m  are  also  of  importance.  Set  =  0  in  Eq.  (12) 
and  integrate,  and  substitute' for  G  and  as  noted  above.  The  height  z  where 
m  is  a  maximum  is  the  solution  of  Eq.  (IB).  . 


el°"4(z0  -  Z) 


30 


60  -  3  X  10 


iT: 


(15) 
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/Solutions  represented  by  Eqs.  (14)  and  (15)  are  plotted  in  Fig.  4.  The 

height  along  each  m-  curve  where  =  0,  the  compensation  level,  is  for 

dz 

each  condensation  level  z  ,  the  height  where  m  =  -G,  i.e.,  where  the 

c  9z 

tendency  of  the  generating  function  to  increase  the  density  of  condensate  is 
exactly  compensated  by  the  decreasing  density  of  ascending  air.  Above  this 
level  of  compensation,  cloud  density  decreases  with  ascent  and  increases 
with  descent. 

The  generating  function  in  this  model  is  zero  at  and  above  10  Km. 

Above  10  Km,  therefore,  the  curves  of  Fig.  4  can  be  extended  by  use  of  the 
relationship 

-4  4 

m(z)  =  m(10Km)  e10  (10  "  z)  ■■  (16) 

which  is  similar  to  that  for  the  other  atmospheric  constituents  not  bound  by  the 
particular  laws  applicable  to  waiter  substance. 

..  This  analysis  may  explain  some  of  the  vagaries  of  the  association  of  high 
cloudiness  with  vertical  ' air  currents.  When  vertical  currents  carry  an  amount 
of  cloud  above  the  level  of  compensation  for  that  cloud  density,  the  cloud  must  be 
very  persistent.  Further  ascent  thins  but  does  not  dissipate  it,  while  descent 
thickens  it.  Ascent  increases  the  size  of  individual  particles  of  condensate,  but 
this  is  more  than  offset  by  their  movement  apart.  The  converse  holds  for  descent 
of  cloud  between  levels  above  the  height  of  compensation.  High  level  horizontal 
convergence  and  descending  motion  of  a  cloud  layer  isolated  at  levels  above  the 
height  of  compensation  by  effects  of  wind  shear,  could  be  associated  with  increas¬ 
ing  cloud  depth  and  dene  tty  instead  of  the  opposite  so  often  assumed,  Of  course, 
Fig.  4  is  inadequate  for  evaluating  the  probable  importance  of  this  effect  in  the 

v>  .  '  " 

real  atmosphere,  since  the  model  generating  function  is  only  a  rough  approxima¬ 
tion  to  the  real  one  at  the  higher  lev  Furthermore,  the  cloud  water  content  is 
in  nature  not  controlled  by  the  generating  function  and  air  density  gradient  alone, 
but  by  these  with  the  depleting  effects  of  precipitation.  Observations  show  that 
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15 


■I 


the  values  of  Fig.  4  are  much  larger  than  usually  occur.  Therefore  the  effect 
of  compressibility,  which  is  proportional  to  m,  must  really  be  of  less  con¬ 
sequence  than  Fig.  4  implies.  During  the  next  quarter,  a  model  generating 
function  which  more  accurately  represents  the  real  tropical  atmosphere, 
especially  at  high  levels  where  G  is  very  small  but  has  great  consequence  for 
cirrus  cloudiness,  will  be  applied  to  an  extension  df  this  study. 


.4  Example  Computer  Solution  of  a  One-dimensional  Case  with  Cloud  Conversion 


and  Accretion 

Quarterly  Progress  Report  No.  1  of  this  contract  discusses  how  the  one¬ 
dimensional  forms  of  Eqs.  (1)  and  (2)  for  an  incompressible  atmosphere  have 
been  programmed  for  solution  by  digital  computer.*  During  the  past  quarter, 
this  program  was  used  to  compute  the  time- dependent  vertical  distributions  of 
precipitation  and  cloud  in  an  initially  saturated  atmosphere.  The  interaction  ■ 
terms  in  Eqs.  (1)  and  (2)  which  are  of  importance  here  are  , 

dM"  dm  .  .  -3  -1  " 

~dtT  or  “  dt"  =  +  ^  m  “  sec  • 


dM  dm  n' r\c  v  1/1*^  0.125  u *  —  - 87 5  —3  .  —1  /10. 

—  or  -  —  =  6.96  X  10  En_  riiM  gm  m  sec  ,  (18) 

at  at  ••..  U  //  . 

a  „  ' 

Eq.  (17)  refers  to  autoconversion  of  cloud  to  precipitation.  Equation  (18)  refers 
to  the  collection  of  cloud  by  precipitation. 

Also,  in  Eq.  (2),  the  fall,  speed  V  is  taken  as  the  fall  speed  VQ  of  the  median . 
volume  diameter  DQ: 

V  =  VQ  =  -  38.G  iIq  M  m  sec  .  (19) 


*Although  the  discussion  of  Section  1.3  of  this  report  demonstrates  that  the 
atmosphere 's  vertical  variation  of  density  may  be  an  important  factor  for  the  cloud 
density  at  high  altitudes,  omission  of  the  air  density  term  is  justified  in  the  machine 
computations  discussed  here  for  two  reasons:  First,  the  effect  of  compressibility, 
proportional  to  cloud  content,  is  small  in  these  cases  because  the  cloud  content  is 
depleted  by  precipitation  to  rather  low  values.  Second,  omission  of  the  compressi¬ 
bility  term  provides  more  direct  understanding  of  the  role  of  the  others.  It  is 
desirable  to  start  with  a  simple  model  and  to  introduce  refinements  progressively. 
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The  use  of  the  cloud  conversion  term  in  Eq.  (17)  is  not  clearly  supported 

or  refuted  by  the  literature,  but  seems  necessary  if  a  greater  number  of  basic 

equations  is  to  be  avoided.  Squires  and  Twomey  [16],  for  example,  suggest 

that  important  variations  in  cloud  droplet  concentrations  and  colloidal  stability 

may  occur  as  a  result  of  variations  of  cloud  nuclei  content.  In  our  initial 

-3 

selection  of  the  value  a  =  0.5  gm  m  ,  we  have  been  guided  by  the  literature 

-3  ' 

that  describes  cloud  water  observations  up  to  about  1  gm  where  precipitation 

was  not  previously  or  subsequently  observed  [1,  15,  17].  The  first  choice  of 

-3  -I 

K  =  10  sec  insures  that  substantial  precipitation  exists  in  the  cloud  20 

*  „  3 

minutes  after  the  cloud  water  content  exceeds  0.5  gm/m  .  The  magnitude  of 

applicable  in  nature  requires  further  clarification  by  reference  to  theoretical 

studies  such  as  [12].  r  -  ...  — 

The  derivation  of  Eq.  (19)  is  discussed  in  Report  No.  1.  This  equation 

defines  the  fall  speed  of  the  median-volume  precipitation  particle  in  the 

exponential  number  distribution  comprising  water  content  M  and  containing 

number  nQ  per  unit  size  range  of  zero-diameter  particles.  It  is  a  little  crude 

to  assume  as  done  here  that  the  vertical  flux  is  given  by  approximating  the 

terminal  fall  speed  of  all  the  precipitation  particles  by  the  fall  speed  of  the 

/median  particle.  This  has  been  shown  by  plotting  the  equation 

„  .  D  ,3  -3.75  D/D.  “  A' 

M  ■»  M^— )  e  P  .  (20) 

which  relates  the  normalized  mass  distribution  of  particles  in  the  Marshall- 
Palmer  distribution  to  size  — .  Slightly  more  than  80%  of  the  water  content  of 
the  exponential  distribution  is°represented  by  drop  diameters  D  such  that 
.4Dq  <  D  <  1  -6Dq,  i.e.,  a  factor  4  in  drop  diameters.  Since  fall  speed  is  propor¬ 
tional  to  the  square  root  pf  drop  diameter,  there  is  a  two  fold  variation  of  fall 
speed  among  the  particle  sizes  which  comprise  80%  of  the  water  content  centered 
about  Dq.  In  other  words,  fall  speeds  of  80%  of  the  water  content  are  in  the 
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Height  IOzm  Height  10  m 


Fig.  5a.  Vertical  distributions  of  precipitation  content  at  1000,  2000,  and  3000 
seconds  in  the  model  defined  by  Table  0,  j 

Fig.  5b.  Vertical  distributions  of  cloud  content  at  1000,  2000,  and  3000  sc<  o".'  s 
in  tho  model  defined  by  Table  0.  1  ! 
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TIME  IN  SECONDS 


HEIGHT 

IN 

METERS 

200 
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1260 

1772 

2000 
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Preclp, 
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Roln 

Starts 

Moximum 

Accretion 
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Cloud 

Woter 

Moximum 

Surfoce 

Preclp. 

Updraft 

Stops 

End  of 
Calculation 

6000 

• 

5700 

m  g/m3  9  Cloud  Water  Content 

g/m3  *  precipi);ation  Content  and  Fall  Speed 
V  m/scc  ▼ 

g/m3  sec  ■  Accretion  of  Cloud  by  Precipitation 

5400 

- 

5100 

4800 

| 

" 

1 

4500 

4200 

■—  ' 

__  . 

-  .  .  .  ■■ 

,, 

•  3.93-I0-1 

3900 

•  761  •  I0*1 

•7.380-I0"' 

3600 

3300 

♦  8. 26  dO-! 

' 

3000 

•  .  _ 

- 

2700 

.. 

2400 

•• 

■  6.60’ 10"* 
*1.35- 10" 2 
▼  3.9 

•  8.10  - 10“ 1 

■  4.00'  I0‘5 

2100 

•2.360,IO'1 

—  ^ - 

1800 

■5.58d0"4 

Al,l6'l0"i 
T  51 

■  I,40'I0'3 

■  1.22'  I0"3 

■1.02- 10“3 

1500 

1200 

■■ 

900 

ii 

..  ••  " 

600 

,, 

> 

*5,50  ■  10-' 
▼  6.2 

il 

300 

1 

11 

0 

’i 

Roln 

Storts 

*6,63' 10-' 
*  6,4 

*6.58' 10-i 
v  6.3 

*1.01  - 10-' 
v  5,0 

Fit*.  G.  Distribution  in  height  and  time  of  significant  space  maxima  in  [the  cloud- 
prcelpitatlon  model  defined  by  Table  G.  Note  that  the  precipitation  maxinjjum  increases 
as  it  descends;  the  cloud  maximum  doorcases  as  it  ascends.  After  the  updraft  ceases 
at  2000  seconds,  the  surface  precipitation  rate  decline^  rapidly  but  the  clbud  aloft  i  . 
relatively  persistent.  r  il 
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The  figures  show,  that  the  cloud  water  content  first  attains  the  value  0.5 

3 

gm/m  required  for  initial  precipitation  development,  in  the  lower  atmosphere 
slightly  above  the  level  where  wG  is  a  maximum.  Once  precipitation  forms,  it 
increases  rapidly  by  collecting  cloud.  At  early  times,  both  cloud  and  precipita¬ 
tion  have  their  maxima  at  the  same  height;  cloud  conversion  and  accretion  pro¬ 
ceed  most  rapidly  at  that  height.  The  cloud  collection  process  is  responsible 
for  the  marked  depletion  of  cloud  in  the  lower  atmosphere,  for  the  steady  down¬ 
ward  shift  of  the  precipitation  maximum  and  the  upward  displacement  of  the 
cloud  maximum.  Note  in  Fig.  6  that  the  precipitation  maximum  progresses 
toward  the  surface  at  a  rate  of  about  2  m/sec  although  the  fall  speed  VQ  is 
typically  6  m/sec.  The  difference  is  important  for  the  interpretation  of  radar 
observations.  The  radar  reflectivity  defined  by  this  model  has  maxima  and 
minima  at  the  same  places  as  M. 

The  solution  of  this  problem  shown  in  Fig.  7  ^hould  be  compared  to  the 
case  shown  as  Fig.  1  in  Report  No.  1  of  this  contract.,  The  parameter  values 
in  the  two  cases  are  identical  except  for  the  distribution  of  G,  a  linear  function 
of  height  in  the  present  model  and  a  constant  in  the  former.  The  average  value 
of  G  is  the  same  in  the  two  cases.  Notice  the  similarity  of  the  curves  in  Fig.  1 
of  Report  No.  1  and  Fig.  7  of  this  report.  It  is  perhaps  surprising  that  the 
solutions  are  not  more  affected  by  such  a  great  change  of  the  distribution  of 
the  generating  function. 

Note  that  the  solutions  are  characterized  by  an  initial  pulse  of  precipita¬ 
tion  and  cloud  content,  then  their  decline  to  lesser  values.  The  pulse  is  due  to 

-3  ''  '  -1 

the  choices  of  a  =  0.5  gm/m  and  =  .001  sec  ,  values  which  with  the 
other  parameters,  lead  to  considerably  higher  cloud  contents  than  can  exist  in 
equilibrium  with  generation  and  accretion  processes.  The  cloud  distribution  „ 
after  precipitation  is  initiated  can  hardly  resemble  the  cloud  distribution  before 
The  onset  of  precipitation  marks  an  alteration  of  the  interplay  of  processes 
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affecting  cloud  distributions,  which  must  persist  until  the  circulation  ceases 
and  all  precipitation  has  fallen  from  the  cloud. 

Note  in  Fig.  6  that  precipitation  spreads  to  levels  above  4,800  meters 
where  the  cloud  water  content  never  reaches  the  value  a  =  0.5.  Such  spread 

of  precipitation  is  an  inherent  feature  of  the  finite  difference  computations. 

\\ 

The  precipitation  content  at  any  height  z  at  time  t  +  At  is  the  mean  of  the 
values  at  z  +  Az  and  z  -  Az  at  time  t,  plus  the  changes  due  to  advection 
and  interaction  processes  computed  for  the  time  interval  At.  This  computational 
method  is  necessary  to  guarantee  that  the  answers  are  both  computationally 
stable  and  better  approximate  the  true  solution  to  the  differential  equations  as 
At  and  Az  become  smaller.*  The  effects  of  averaging  appear  partly  as  a  very 
slow  depletion  of  cloud  by  precipitation  even  above  levels  where  the  spontaneous 
conversion  of  cloud  to  precipitation  is  active  and  even  when  VQ  is  markedly 
greater  than  maximum  updrafts.  This  manifestation  probably  has  a  natural 
analogue,  since  in  any  precipitation,  however  intense,  there  must  exist  some 
drops  whose  fall  speeds  are  less  than  the  updrafts,  and  which  are  carried  upward 
with  the  updrafts. 

1.5  Relationship  of  Precipitation  Rate  to  the  Updraft  Intensity  and  to  the  Rate 
of  Cloud  Conversion  , 

The  parameter  values  in  Table  7'have  been  used  to  test  the  effect  of 

variations  of  updraft  strength  and  rate  of  cloud  eonversion  on  the  distributions 

of  cloud  and  precipitation. 

♦Tests  have  shown  how  the  solutions  vary  with  choice  of  Az  and  At.  In  the  n 
problem  discussed  here,  a  typical  cloud  or  precipitation  value  at  a  given  time  " 
is  shifted  3%  for  a  change  of  space  discretization  from  300  to  150  meters,  3% 
more  for  a  change  from  15p  to  60  meters,  but  just  1%  additional  for  a  shift  of 
Az  from  60  to  30  meters  (200  steps).  The  use  of  100  space  steps  seems  a 
suitable  compromise  between  accuracy  and  economy. 


TABLE  7* 

PARAMETER  VALUES  FOR  STUDY  OF  EFFECTS 
OF  UPDRAFT  STRENGTH  AND  RATE  OF  CLOUD  CONVERSION 
ON  VERTICAL  DISTRIBUTION  OF  PRECIPITATION  AND  CLOUD. 

Case  a)  w  =  0.5  m/sec  K,  =  10  ^  see  * 

max  1  -3  -i 

Case  b)  w  =  0.5  m/sec  K„  =  10  Sec  1 

..  max  1  '  „5  _1: 

Case  c)  w  =  0.5  m/sec  KL  =  10  sec 

-  ;  max  1  _3  ^ 

Cases  d  and  e)  w  =  2.0  m/sec  K..  .=  10  sec  1 

'  max  1 

^Parameters  not  identified  here  have  been  assigned  the  values  given  in 
Table  6  in  Section  1.4  of  this  report.  - 

The  surface  precipitation  rates  given  by  these  models  are  depicted  in 

Figs.-  8  and  9.  Figure  8  shows  that  the  onset  time  of  precipitation  increases  .. 

about  10  minutes  for  a  decrease  of  by  four  orders  of  magnitude;  the  peak 

rate  of  rainfall  increases  only  slightly  for  this  enormous  change  of  K^.  This 

type  of  finding  may  have  some  implications  for  cloud  seeding  efforts  but 

should  not  be  considered  definitive  at  this  stage  of  the  work.  The  result  seems 

"  "5 

astonishing  when  one  considers  Eq.  (17)  alone.  With  =  10  ,  about  .01  of 
the  cloud  water  in  excess  of  the  amount  a  spontaneously  becomes  precipitation 
after  1000  secs.f  yet  appreciable  precipitation  at  the  ground  occurs  in  this 
model  only  600  seconds  later  than  in  the  case  where  spontaneous  conversion 
of  cloud  water  is  10,000  times  faster.  This  result  testifies  to  the  effectiveness 
of  accretion,  Eq.  (18),  even  when  the  amount  of  precipitation  is  very  small. 

This  result  should  not  be  considered  indicative  of  unimportance  of  the 
rate  of  cloud  conversion  in  natural  processes.  Since  a  nonlinear  relation  between 
cloud  conversion  and  collection  governs  the  early  rate  of  precipitation  increase, 
substantial  changes  are  probably  associated  with  variations  of  when  E  and 

fThe  amount  of  cloud  water  in  this  model  which  becomes  precipitation  in  1000 
seconds  by  virtue  of  spontaneous  cloud  conversion  alone  is  about  .005  gm/m3. 
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Fig.  8.  Time  dependence  of  surface  precipitation  rate  as  for  three  valaes  of  K,  the  rate  at  which  cloud  content  in 
excess  of  the  amount  a  changes  spontaneously  to  precipitation.  The  sluggish  variation  of  surface  precipitation  with 
very  large  change s  of  K  indicates  the  great  importance  of  the  accretion  process  to  comprehensive  evaluation  of  the 
nex  effects  of  cloud -precipitation  interactions.  (Case  b  computations  terminated  at  2600  secs.) 


Rate  of  Rainfa  H  i  mm  /  hr. 


iIq  lie  within  particular  ranges.  The  Eqs.  (17)  and  (18)  can  be  combined  to  form 
a  pair  of  differential  equations  with  a  straight-forward  analytic  solution  which 
clearly  illuminates  the  role  of  the  equations'  coefficients  in  precipitation  initia¬ 
tion.  These  solutions  are  now  being  examined  and  will  be  discussed  in  the  next 
report.  Other  theoretical  work  in  precipitation  physics  [13,  for.example] , 
is  also  being  studied  in  order  to  derive  estimates  of  K^,  n^,  and  E  which  best 
model  natural  conditions.  The  variation  of  collection  efficiency  with  size  of 
precipitation  particles  affects  both  the  constant  coefficient  and  the  exponent  of 
M  In  Eq.  (18). 

At  least  one  "other  property  of  the  solutions  shown  in  Fig.  8  is  worthy  of 
note.  At  4000  seconds,  the  precipitation  rates  at  the  ground  (and  the  distribution 
of  cloud  and  precipitation  aloft)  are  almost  exactly  the  same  in  the -three  cases. 
However,  the  precipitation  accumulated  by  4000  seconds  is  11.96  mm  when  ., 

-1  -1  _5  _,i 

=  10  sec  and  declines  to  just  10.08  mm  when  ■  10  sec  .  The 
smaller  value  of  is  associated  with  longer  initial  duration  of  large  amounts 
of  cloud  in  the  upper  atmosphere  where  it  is  removed  by  horizontal  divergence 
of  the  wind.  \ 

Figure  9  shows  the  precipitation  rate  at  the  ground  in  three  cases  of 
Table  7.  Case  b  in  Fig.  9  is  the  same  as  Case  b  in  Fig.  8.  The  value  of 
in  Case  d  is  the  same  as  in  b,  but  the  maximum  updrafts  in  Case  d  are  2m/sec 
and  stop  after  600  sec.  Thus  the  amount  of  condensate  produced  in  2000  seconds 
in  Case  b  equals  that  produced  In  600  seconds  in  Case  d.  This  amount  on  a  level 
surface  would  have  a  depth  of  8  mm.  Now  consider  the  precipitation  totals 
received  at  the  ground  In  two  cases.  By  2000  seconds,  Case  d  produces  6.3  mm 
of  precipitation  at  the  surface;  about  .2  mm  more  falls  after  2000  seconds  for  a 
total  of  6.5  mm.  Case  b,  on  the  other  hand,  produces  oii'ly  3.5  mm  of  precipita¬ 
tion  by  2000  seconds,  but  2.5  mm  falls  after  that  time  when  the  updrafts  stop, 

3  for  a  total  of  6.0  mm.  Neither  case  produces  as  much  as  is  condensed,  since 
some  condensate  still  remains  as  cloud  aloft  when  computations  are  terminated,!! 
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and  some  condensate  is  lost  by  divergence  at  high  levels  during  the  time  period 
before  precipitation  starts. 

The  greater  efficiency  of  Case  b  than  Case  d  can  be  understood  in  terms 
of  the  discussion  given  in  connection  with  Eq.  (7),  Section  1.1.  First,  the  low 
level  convergence  which  acts  to  increase  precipitation  at  updraft  cores  at  the 
expense  of  amounts  away  from  the  core,  is  effective  only  if  the  precipitation 
descends  through  updrafts.  In  Case  d,  the  updrafts  stop  before  appreciable 

it  •  ..  ■  .. 

precipitation  reaches  the  ground.  Second,  in  the  strong  updraft  case,  cloud  is 
lifted  higher  before  it  changes  to  precipitation;  therefore  a  greater  amount  of 

cloud  is  removed  from  the  updraft  core  by  divergence  in  Case  d  than  in  Case  b. 

ii 

To  the  extent  that  cloud  is  carried  to  greater  heights  by  stronger  updrafts  and 

more  spread  out  there  by  accompanying  .strong  horizontal  divergence,  instead 

of  descending  to  the  ground  as  precipitation  in  the  Saturated  updrafts ,  the 

stronger  updraft  cases  are  less  efficient  producers  of  precipitation  than  are  the 

weak  updraft  models.  Accurate  evaluation  of  these  effects  will  be  possible  when 

two-  and  three-dimensional  models  are  programmed  for  the  computer. 

In  Case  e,  2  m/sec  updrafts  continue  to  3000  secs.  Note  that  the  steady 

precipitation  rate  in  this  case  is  somewhat  more  than  four  times  greater  than 

the  precipitation  rate  in  Case  b.  That  the  steady  precipitation  rate  is  not  exactly 

proportional  to  the  updraft  intensity  is  related  to  variations  with  updraft  strength 

"of  the  shape  of  the  vertical  distributions  of  cloud  and  precipitation  and  to  changes 

of  the  proportional  contribution  of  the  divergence  term  in  Eq.  (7). 

It  is  particularly  worthy  of  note  in  these  solutiohs,  as  in  others  not  repro- 
»  .  '  -  .  "  "  ••  •• 
duced  here,  that  the  transient  effects  manifested  in  the  surface  precipitation 

rates  by  the  simulated  cloud  physics  processes  are  rather  small.  The  surface 

precipitation  rate  is  mightily  dependent  on  the  updrafts,  however.  Part  of  our 

'  task  during  the  next  interval  will  be  to  learn  if  there  are  parameter  combinations 

“  -  K\ 

which  produce  larger  oscillations  or  other  transient  effects  with  steady  updrafts 
than  the  weak  variations  computed  thus  far. 
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1.6  Comment  on  Mathematical  Background 

The  mathematical  background  of  the  equations  employed  in  this  study 
has  been  investigated  starting  with  the  one-dimensional  cloud-free  model 
where  the  fall  speed  of  condensate  is  constant  [7,  8].  The  existence  and 
uniqueness  of  the  solutions  presented  and  discussed  in  the  references  from 
a  physical  standpoint  is  rigorously  demonstrated  below. 

The  one-dimensional  cloud-free  equation  is 

.  IT  +  <v  +  *)  if  “  *G 

in  which  V  is  constant; 

w  -  az(H  -  z) 

and  G  is  a  differentiable  function  of  z  alone.  .  Equation  (21)  is  to  be  solved  in 


the  strip 


S:  0  S  u  S  H; 


0  s  t  <  +  • 


subject  to  the  conditions 

M<a'0)“l  (28) 

M(H,  t)  -  0  ' 

and 

V  +  w  <  0  ,  „  (24) 

The  physical  meaning!  of  Eq.  (24)  is  that  updrafts  are  everywhere  less  than  the  < 
precipitation  falling  t  peed.  The  problem  as  posed  has  a  continuous  unique 
solution,  proved  as  f<  llCWs.  The  characteristic  equations  of  (1)  are 

•  S-  u  g-V.W,  ® -WO  (25) 

Through  each  point  of  the  interval 

t  =  0,’  OsjSH;  M»0,  (26) 

will  pass  a  unique  continuous  solution  of  Eq.  (25)  subject  to  the  initial  conditions 


t  =  M  =  0 


(0  s  zQ  s  H) 


9" 
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and  through  each  point  of  the  half-line 

z  =  H,  0  s  t  <  +  *>,  M  =  0 

will  pass  a  unique  solution  subject  to  the  initial  conditions 
z  =  H,  M  =  0 

t  =  tQ  -  "  (°  s  tfl  <  +  «). 

The  uniqueness  and  existence  of  the  solutions  follow  from  standard  existence 
theorms  [5]  since  the  right  sides  of  Eq.  (25)  are  continuous  and  Lipschitzian. 

In  the  t  -  z  plane,  the  characteristic  base  curves  for  Eq.  (21)  are 

(30) 

These  base  curves  cover  the  strip  S 
in  the  sense  that  through  each  point  in  the  strip  passes  one  and  Only  one  base 
curve.  Each  of  these  curves  intersects  one  or  the  other,  but  not  both,  of  the 
intervals  (26)  and  (28)  in  one  and  only  one  point.  These  facts  follow  from 
Eq.  (24)  and  the  Lipschitzian  character  of  V  +  w.  With  each  point  on  a  base 
curve  in  the  strip  S  there  is  associated  a  unique  value  of  M  determined  from 
(25)  with  the  proper  initial  point.  The  M  associated  in  this  manner  with  each 
point  in  the  strip  S  determines  a  unique  continuous  solution  of  (21)  which 
satisfies  (23)  as  follows  from  the  standard  theory  of  linear  partial  differential 
equations  in  two  independent  variables  [3]. 

The,  method  of  solution  is  summarized  a#  follows.  To  determine  the 
value  of  M  at  the  'point  (z^  t^)  in  the  strip  S,  solve  (30)  to  determine  the 
characteristic  base  curve  which  goes  through  (z^,  t^).  Call  this  curve  C^. 

Find  the  point  of  intersection  of  Cj  with  either  (26)  or  (28)  , as  the  case  may  be. 
Call  t)iis  point  (Zq,  t  ).'  Determine  M  from'  (25)  by  solving  for  M  as  a  function 
of  z  from  • 

dM  wG  ” 

dz  “  V  f  w  <31) 

with  M(z^)  =  0,  and  setting  z  =  z^.  H  - 


determined  from 

dz  „  ,  - 
•j-V.w  ■ 

together  with  proper  initial  conditions. 


(28) 

(29) 
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One  should  notice  that  the  characteristic  base  curve  which  goes  through 
the  point  t  =  0,  t  =  H  divides  the  strip  S  into  two  regions:  one  to  the  right 
(t  increasing)  and  one  to  the  left  (t  decreasing)  of  the  base  curve.  For  points 
in  the  strip  which  lie  to  the  right  of  this  curve,  the  value  of  z0  of  the  above 
paragraph  is  always  H  and  the  solution  of  (31)  is  independent  of  t.  That  is, 
the  solution  of  (21)  under  conditions  (23)  and  (24)  in  the  strip  S  attains  &  steady 
state  for  which  5M/8t  =  0  for  z  fixed. 

An  interesting  variation  of  the  problem  considered  above  iB  defined  by 
conditions 

V  +  w  <  0  0  s  a  <  h^  hg  <  z  s  H 

V  +  w  =  0  z  .=  hj,  hg  .  (32) 

V  +  w  >  0  hj  <  2  <  h2 

The  conditions  (32)  mean  that  the  maximum  updraft  exceeds  the  precipitation 
fall  speed.  Application  of  the  previous  reasoning  to  each  of  the  three  strips 
in  the  t  -  z  plane  determined  by  (32)  for  t  as  0  leads  to  the  conclusion  that 
there  exists  a  unique  continuous  solution  of  (21)  for  0  s  z  s  H,  t  a  0,  which 
satisfies  (23)  under  coiidltions  (32). .  Consideration  of  the  characteristic  base 
curves  of  (21)  for  this  problem  shows  that  for  any  fixed  z  >  h2  there  is  a 
time  beyond  which  the  solution  is  steady.  For  any  z  >  h^,  there  is  no  steady- 
state  solution. '  this  case  is  discussed  from  a  physical  point  of  view  in  [8] . 

The  question  of  existence  and  uniqueness  for  the  quasi-linear  hyperbolic 
system  describing  clouds  and  rain  cannot  be  answered  at  present.  It  is  a 
reasonable  conjecture,  that  a  unique  continuously  differentiable  solution  does  , 
not  exist  for  conditions  of  M,  m  »  0  initially.  R.  Courant  [2]  has  recently 
established  the  existence  of  a  continuously  differentiable  solution  for  the 
hyperbolic  system,  +  A(x,  t,  u)ux  +  B(x,t,  u)  «=  0,.  with  initial  data 
u(x,  0)  =  #(x)  on  some  interval  of  the  x-  axis.  His  method  of  proof  requires 
that  A(x,t,u)  ahd  B(x,  t,  u)  be  continuously  differentiable  in  the  domain  of 


31 


x,  t,  u  under  consideration.  Our  problem  fails  to  satisfy  this  condition. 

Since  Courant's  proof  follows  from  the  weakest  hypotheses  deemed 
reasonable  for  success  in  solving  the  problem,  it  is  felt  that  a  unique 
solution  for  the  equations  of  clouds  and  rain  must  perforce  be  a  generalized 
solution  which  allows  discontinuities  in  the  derivatives  of  the  solution  veotor. 

We  believe  that  a  solution  of  this  type  might  be  established  along  the  lines 
followed  by  Kuznecov  and  Rozdestvenskii  [10] .  However,  due  to  the  rather 
abstract  and  purely  mathematical  nature  of  this  problem  it  is  felt  that  no  . 
further  effort  should  be  devoted  to  it  on  this  contract. 

The  incomplete  theory  of  the  techniques  used  in  this  investigation  is  not 
satisfying  to  a  mathematician.  However,  this  fact  is  not  cause  for  rejection  of 
the  results  of  the  computations.  The  physical  reasoning  underlying  the 
mathematicai  formuiation  of  this  problem  and  the  numerical  results  of  these 
and  similar  calculations  of  others  (e.g.,  P.  D.  Lax  [3])  imply  strong  credibility 
for  the  results.  The  semi-empirical  approach  used  here  is  typical  in  the 
development  of  new  methods.  As  R.  D.  Richtmyer  [4]  has  said,  "...  if  we  were 
to  wait  for  convergence  proofs  and  error  estimates  for  the  new  methods,  most 
of  the  computers  now  in  use  in  technology  and  industry  would  come  grinding  to. 
a  halt." 

1.7  The  Difference  Between  the  Horizontal  Speed  of  Raindrops  and  the 

Horizontal  Wind  Speed. 

The  precipitation-cloud  models  discussed  in  this  report  are  based  in 
part  on  the  assumption  that  the  horizontal  speed  of  condensate  is  the  same  as  the 
horizontal  speed  of  the  air.  The  validity  of  this  assumption  must  decrease  with 
increasing  size  of  precipitation  particles;  a  fair  test  of  this  assumption  is  there¬ 
fore  given  by  examination  of  an  equation  of  motion  for  a  large  raindrop. 

The  horizontal  resisting  force  on  a  spherical  raindrop  is  given  by 

CDR 

f ,  =  6  rpr(u  -  v)~ j-*  (83)  „ 
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where  ft  is  the  viscosity  of  air,  r  is  the  radius  of  the  drop,  u  is  the  horizontal 
velocity  of  the  air,  v  is  the  horizontal  velocity  of  the  drop,  C^  is  the  drag 
coefficient,  and  Rg  is  the  Reynold's  number  [6,  p.  121].  Consider  a  drop 
falling  at  terminal  speed  through  a  layer  of  vertical  shear  of  the  horizontal 
wind.  The  drop  is  accelerating  horizontally  due  to  the  force  associated  with 
a  difference  between  its  horizontal  velocity  and  the  horizontal  air  motion, 
This  force  given  by  Eq.  (33)  is  equal  to  the  mass  of  the  drop  times  its 
horizontal  acceleration,  i.e., 


a  ..  °DRe  4ir  3  dv 

Qirfir  (u  -  v)  =  -j  r  P  dt’  (34) 

where  p  is  the  density  of  water.  Consider  a  drop  of  radius  0.24  cm.  The 

-4  -1  -1 

quantity  C^R^/24  «  77  and  fi  ~  1.8  x  10  gm  cm  sec  .  Substitute  these 
values  in  (34)  to  obtain 

=  1.09  (u  -  v).  (35) 


Nothing  of  importance  to  this  development  is  lost  by  considering  the  case  where 
the  shear  is  uniform  and  the  difference  between  the  horizontal  wind  speed  and 
drop  speed  is  constant.  The  drop's  falling  speed,  about  907  cm  sec  *  and  the 
shear  8u/9z  *  B  determines  du/dt,  the  time  rate  of  change  of  u  (and  v)  at  the 
drop.  „ 


du  _  „dv  dz  Bu 
dt  dt  ~  dt  8z 


907B. 


(36) 


Substitute  for  dv/dt  in  (35)  from  (36)  to  get 


u  -  v  = 


907B 

1.09 


(37) 


When  B  is  0.02  sec  (20  m/sec  km,  a  very  large  value),  u  -  v  is  only  16.64 

j 

cm/sec,  less  than  1%  of  representative  horizontal  wind  speeds, 


So  long  as.  u  -  v  is  small  compared  with  the, terminal  fall  speeds,  the 
method  employed  here  yields  an  accurate  measure  of  u  -  v.  The  generalization 
to  variable  shear  and  variable  u  -  v  is  obvious  but  of  little  practical  interest 
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to  the  work  of  the  contract  since  u  -  v  is  so  small.  The  work  of  Das  [4, 

Table  1  and  app,] ,  brought  to  our  attention  after  that  above  had  been  com¬ 
pleted,  leads  by  a  different  path  to  the  same  conclusion. 
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2.0  CONCLUSION; 

The  basic  formulation,  Eqs.  (1)  and  (2),  IS  a  valuable  tool  for  Investiga¬ 
tion  of  the  macroscopic  effects  of  cloud  physics  processes  and  of  fundamental 
properties  of  associations  between  wind  and  condensate  distributions.  Solu¬ 
tions  indicate  much  about  the  extent  of  control  exercised  by  the  cloud  physics 
processes  on  the  shapes  of  cloud  and  precipitation,  distributions  and  oil 
precipitation  intensity , 


/ ' 
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3.0  PROGRAM  FOR  NEXT  INTERVAL 


The  investigative  program  during  the  next  interval  should  include  study 
in  the  following  areas,  not  all  of  which  will  be  concluded. 

1.  Accurate  determination  of  the  profile  of  saturation  water  vapor 

density  in  the  tropical  atmosphere,  especially  at  high  levels,  and  examination 
of  the  implications  of  this  distribution  for  the  persistence  of  high  cloud 
(Section  1.3).  ■■  »' 

2.  Study  of  literature  pertinent  to  selection  of  cloud  conversion  co¬ 
efficients  and  collection  efficiencies  representative  of  the  real  atmosphere 
(Sections  1.4  and  1.5). 

3.  Examination  of  the  simultaneous  analytic  solutions  of  equations 
modeling  the  time  dependence  of  m  and  M  in  the  presence  of  cloud  conversion 
and  accretion  processes,  and  analysis  of  the  role  of  cloud  conversion  and 
accretion  coefficients  in  shaping  precipitation  development  (Section  1.6). 

4.  Evaluation  of  the  roles  of  the  shape  of  the  distribution  of  precipitation 
particles  and  of  the  critical  cloud  water  content  for  the  intensity  and  time 
dependence  of  surface  precipitation  and  vertical  profiles  of  precipitation  and 
cloud.  This  work  will  proceed  by  examination  of  the  computer  solutions  for 
selected  combinations  of  nQ  and  a  (Section  1.6). 

5.  Preparation  of  a  computer  program  for  two  dimensional  and  radially 

symmetric  models  with  provisions  for  computation  of  cloud,  precipitation,  and 
vapor  budgets  within  the  atmospheric  volume.  „ 
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4.0  IDENTIFICATION  OF  PERSONNEL 
4.1  Extent  of  Participation 


Name 

title 

•Via  .  •  _ 

Total  hours  worked 
during  the  quarter 
(approximate) 

E.  Kessler 

Principal  Investigator 

100 

P.  J.  Feteris 

Research  Scientist  (Meteorologist) 

330 

E.  A.  Newburg 

Research  Scientist  (Mathematician) 

136 

C.  P.  Robson 

Research  Associate  (Programmer) 

30  \  \ 

Q.  Wiokham 

Research  Associate  (Programmer) 

130 

P.  Duohcrw 

Senior  Research  Aide 

70 

Secretarial,  administrative,  and  drafting  assistance  was  also  provided. 

Mr.  Feterls  joined  the  project  on  August  0,  1962. 

4.2  Biography  erf  Plater  J.  Feterls 

Mr.  Feteris  is  a  research  Scientist  in  the  Atmospheric  Physics  Division 
of  The  Travelers  Research  Center,  Ino.  He  is  taking  part  in  the  analysis  and 
development  of  cloud  and  precipitation  models  and  is  also  engaged  in  a 
statistical  study  to  relate  the  statistical  properties  of  radar  patterns  to  pattern 
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From  1967  to  1962  he  was  Junior  Scientific  Officer  at  the  Royal  Netherlands 
Meteorological  Institute,  where  he  made  weather  analyses  and  issued  shipping 
forecasts.  , 
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and  he  received  his  MS  in  Meteorology  from  The  University  of  London  in  I960. 
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